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ABSTRACT 

Binary  sequences  have  had  application  in  communication 
systems  for  many  years.   Shift  registers  have  been  used  in 
their  generation,  because  of  the  ease  and  economy  of  their 
operation.   For  certain  applications,  nonlinear  feedback 
functions  are  used  by  shift  registers  of  span  n  to  generate 
sequences  of  lengths  up  to  2  . 

The  sequences  of  maximum  length  2   and  their  generation 
are  the  subject  of  this  thesis.   In  particular  the  ways 
of  generating  these  sequences  using  nonlinear  feedback 
shift  registers  and  their  correlation  to  linear  feedback 
shift  registers  are  described.   Complexity  is  the  term 
given  to  the  length  of  the  shortest  linear  feedback  shift 
register  generating  a  maximum  length  2   sequence. 

Games  and  Chan  [Ref.  I]  have  given  considerable  study 
to  the  subject  of  complexity.   Some  of  the  problems  they 
left  are  discussed  further  in  this  paper.   It  will  be  shown 
that  the  complexity  of  a  de  Bruijn  sequence  (S)  is  the  same 
as  the  complexity  of  its  reverse  (r  S)  ,  complement  (5~)  ,  and 
its  reverse  complement  (r  S)  .   Sequences  (S)  for  which  r  S  = 
S~  are  termed  RC  sequences.   It  is  shown  that  RC  sequences 
exist  for  every  odd  n>_3 .   In  addition  a  lower  bound  will  be 
established  for  the  number  of  RC  sequences  occurring  for 
each  odd  n. 
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I.   INTRODUCTION 

A.  APPLICATION 

de  Bruijn  sequences  have  been  the  subject  of  recurring 
interest  since  their  discovery  over  a  century  ago.   Recent 
developments  in  radar  ranging  [Ref.  2:  Ch .  13],  error  correct 
ing  codes,  secure  or  limited  access  code  generators  [Ref.  3: 
pp.  12-14],  and  mathematical  modeling  have  made  application 
of  de  Bruijn  sequences  because  of  their  properties.   The 
binary  shift  registers  readily  model  digital  computers  with 
electronic  states  where  "1"  indicates  the  on  state  and  "0" 
indicates  the  off  state. 

B.  DEFINITION  OF  de  BRUIJN  SEQUENCES 

For  a  positive  integer  n,  a  de  Bruijn  sequence  of  span 

n  is  a  complete  binary  cycle  of  length  2   which  is  a  sequence 

S  =  {s,,  s0,  •••,  s   }  taken  in  circular  order  (s-,  follows 
1    L  7n  l 

s   )  such  that  all  possible  n-tuples  occur  exactly  once. 

2n 
[Ref.  4:  p.  120].   Complete  binary  cycles  are  often  called 

de  Bruijn  sequences  after  the  Dutch  mathematician,  N.  de 

Bruijn,  who  in  1946  proved  the  existence  of  complete  binary 

2n~^-n 
cycles  having  length  2   and  that  they  numbered  2 

de  Bruijn  was  preceeded  by  Flye  Saint-Marie  a  half  century 

earlier  in  an  obscure  publication  [Ref.  5].   As  an  example 

,3-1. 3 

for  n  =  3,  Table  1.1  shows  both  (2      )  de  Bruijn  sequences 
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3 

of  length  eight  (2  ) .   Both  of  these  sequences  are  examined 

to  confirm  that  each  of  the  eight  3-tuples  occur  exactly 
once.   Notice  that  2  elements  of  the  sequence  are  repeated 
as  signified  by  the  bar  (00")  ,  which  was  necessary  on  each 
sequence  to  complete  the  last  two  3-tuples. 

TABLE  1.1 
3-TUPLE  COMPOSITION  OF  de  BRUIJN  SEQUENCES  OF  SPAN  3 


0001110100  0001011100 

000  000 

001  001 

011  010 

111  101 

110  011 

101  111 

010  110 

100  100 


There  are  2n  different  possible  starting  points  for 
each  de  Bruijn  sequence;  however,  each  cyclic  permutation 
is  considered  equivalent.   It  is  often  a  matter  of  conveni- 
ence to  start  each  sequence  in  a  canonical  way  with  the  all 
zero  or  all  one  n-tuple.   This  reduces  confusion  when  working 
with  a  number  of  different  de  Bruijn  sequences. 

C.   PSEUDO- RANDOMNESS 

Since  each  binary  n-tuple  occurs  exactly  once  on  the 
cycle,  the  sequence  models  a  uniform  distribution  with  the 
probability  of  given  n-tuple  occurring  to  be  exactly  /y11- 
The  sequence  is  not  entirely  random,  however,  since  each 
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n-tuple  (a-.  a~  •••  a  )  has  but  two  possible  successors 

(0  a,  a,  •••  a   ,  or  1  a,  a~  • • •  a  ,) .   de  Bruiin  sequences 
■12      n-1       12      n-1  J  n 

satisfy  the  randomness  properties  of  equidistribution  of 
O's  and  l's  and  run  lengths  [Ref.  3:  p.  10]  which  would  be 
expected  from  the  tossing  of  a  fair  coin.   Thus,  de  Bruijn 
sequences  have  a  pseudorandom  property. 

D.  K-ARY  SEQUENCES 

In  general  it  is  possible  to  work  with  k  possible  states 

for  each  position  of  an  n-tuple.   Five  Saint-Marie  in  1894, 

showed  the  existence  and  determined  the  number  of  complete 

,n-l  T,n-1_ 
cycles  for  n-tuples  of  k  characters  to  be  C(k-l)!]    k 

[Ref.  5].   This  report  is  restricted  to  the  binary  case 
(k  =  2)  due  to  its  principal  application  in  electronic  com- 
munications and  computers. 

E.  DEFINITIONS 

For  S  =  {s-.  s2  •••  s-,}    where  s-  e  {0,1}  for  l^i£k  the 
following  are  defined: 

1.  Let  W(S)  denote  the  weight  of  S. 

k 
W(S)  =  Z    s- 
i  =  l  x 

2.  Let  P(S)  denote  the  parity  of  S. 

k 
P(S)  =  Z  s.    (mod  2) 
i  =  l  x 
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3.   Let  1(S)  denote  the  length  of  S. 

l(s)  =  k  where  k  is  the  number  of  positions  in  S 
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II.   SEQUENCE  GENERATORS  AND  OPERATORS 

A.   SHIFT  REGISTERS 

In  practice  shift  registers  are  used  to  generate  de  Bruijn 
sequences.   An  n-stage  shift  register  has  n  memory  registers 
(x, ,  x?  ,  •••,  x  )  which  shift  their  contents  (0  or  1)  to  the 
next  register  upon  command.   For  example,  in  Fig.  2.1,  the 
contents  of  x,  will  transfer  to  x^ ,    x2  to  x,,  etc.,  at  the 
appropriate  time  with  X   serving  as  the  output.   However, 
the  contents  of  an  n-stage  shift  register  would  empty  in 
n  shifts  if  no  input  was  provided  to  x, . 


i 1 


-*i 


i-* 


i i 


OUT 


Fig.  2.1    n-Stage  Shift  Register. 

B.   FEEDBACK  SHIFT  REGISTERS 

Feedback  networks  are  added  to  provide  an  input  to  the 
1st  stage  of  the  shift  register.   Note  that  in  practice 
the  output  can  be  taken  from  any  stage  of  the  shift  register 
In  fact,  the  contents  of  any  stage  is  the  same  as  the  con- 
tents  of  any  stage  is  the  same  as  the  contents  of  the  n 
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stage,  merely  shifted  by  a  certain  amount.   In  Fig.  2.2  the 
feedback  function  f  is  added  to  the  n-stage  shift  register 
to  generate  non-trivial  sequences.   The  coefficients 


aQ,  ax, 


n 


are  0  or  1  with  an  =  a_  =  1 


0 


n 


Y 

4, 

x2 

1 

Xl 

W 

i 

Feedback 
Loop 


OUT 


Fig.  2.2    n-Stage  Feedback  Shift  Register. 


1 .   Linear  Feedback  Functions 

A  feedback  function  f  is  classified  as  linear  if 

the  function  f  (x,  ,  x-> ,  •••.  x  )  is  restricted  to  be  of  the 

1    L  n 

n 

form  f  =  Z  a-  x-  where  the  addition  is  modulo  2  addition. 

i  =  l  1      x 

The  addition  (mod  2)  table  is  given  in  Table  2.1  below. 

Note  that  subtraction  is  equivalent  to  addition  when  opera 

ting  modulo  2. 
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TABLE  2.1 
ADDITION  (MODULO  2) 


® 

0 

1 

0 

0 

1 

1 

1 

0 

2 .   Nonlinear  Feedback  Functions 

If  the  feedback  function  f (x, ,  •••.  x  )  utilizes 

1        n 

multiplication  (mod  2)  and  addition  (mod  2)  then  the  feed- 
back function  is  nonlinear.   (We  also  say  that  the  linear 
feedback  functions  are  vacuously  nonlinear  feedback  functions.) 
Modulo  2  multiplication  is  identical  to  multiplication  in 
the  integers  with  the  restriction  of  the  alphabet  used 
being  the  set  {0,  1}  as  described  in  Table  2.2. 

TABLE  2.2 
MULTIPLICATION  (MODULO  2) 


• 

0 

1 

0 

0 

0 

1 

0 

1 

3 .   Comparison  of  Linear  and  Nonlinear  Feedback  Functions 
Considerable  research  has  been  devoted  to  the  study 
of  linear  feedback  functions  because  of  their  ease  of  analy- 
sis.  Nonlinear  feedback  functions  are  much  more  complicated. 
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However,  the  class  of  sequences  that  can  be  generated  by 
nonlinear  functions  is  much  greater  than  for  linear  functions. 
The  number  of  binary  maximum  length  (2  -1)  linear 

sequences  that  can  be  generated  by  a  n-stage  linear  feed- 

2n 
back  shift  register  (LFSR)  is  approximately  —  [Ref.  7]. 

The  reason  that  the  maximum  length  linear  sequences  only 

have  length  2n-l  is  that  the  zero  state  is  a  fixed  point 

under  the  linear  function.   This  poses  no  problem  since  the 

maximal  length  linear  sequences  are  readily  made  into 

de  Bruijn  sequences  having  length  2   simply  by  adding  a 

zero  to  the  (n-1) -tuple  of  zeros.   The  feedback  function 

involved  then  becomes  a  nonlinear  function. 

By  comparison  the  number  of  binary  nonlinear  full 

length  sequences  (2  )  that  can  be  generated  by  a  nonlinear 

feedback  shift  register  (NFSR)  is  exactly  2     n.   Below, 

Table  2.3  compares  these  numbers  for  3<n<8 .   The  number  of 

maximum  length  sequences  will  be  of  interest  later  in 

Section  III,  where  the  concept  of  complexity  in  generating 

these  sequences  will  be  developed.   Various  algorithms  are 

listed  in  a  survey  by  Fredricksen  [Ref.  5]  for  generating 

de  Bruijn  sequences. 

C.   SUCCESSOR  AND  PREDECESSOR  STATES 

The  contents  of  an  n-stage  shift  register  at  a  specific 
time  is  usually  referred  to  as  its  "state".   As  mentioned 
in  Section  I  when  discussing  the  pseudorandomness  of  de  Bruijn 
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Linear 

Nonlinear 

2 

2 

4 

16 

6 

2048 

6 

67,108,864 

18 

17 
1.4  x  1017 

30 

2.7  x  10°6 

TABLE  2.5 
NUMBER  OF  MAXIMUM  LENGTH  SEQUENCES 

n 
3 

4 
5 
6 
7 
8 

sequences  each  state  has  two  possible  successor  (conjugate) 
states  depending  on  whether  the  feedback  function  generates 
a  "0"  or  "1"  as  the  input  to  register  x, .   Likewise  each 
state  has  two  predecessor  (companion)  states  depending  on 
whether  the  previous  output  was  a  "0"  or  "1".   The  adjacency 
quadruple  in  Fig.  2.3  centers  around  then  (n-1) -tuple  X2#**x 
That  is,  each  input  state  has  two  possible  successors  depend- 
ing on  whether  the  feedback  is  0  or  1,  and  each  output  state 
has  two  possible  predecessors  depending  on  whether  the  pre- 
vious  output  from  the  n   stage  was  0  or  1 . 

In  Table  2.4  we  show  an  example  for  n  =  3  of  the  possible 
successors  and  predecessors  for  each  state. 
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Fig.  2.3   Adjacency  Quadruple. 


TABLE  2.4 
CONJUGATE  AND  COMPANION  STATES  (n=3) 


Predecessors 

State 

Successors 

000-001 

000 

000-100 

010-011 

001 

000-100 

100-101 

010 

001-101 

110-111 

011 

001-101 

000-001 

100 

010-110 

010-011 

101 

010-110 

100-101 

110 

011-111 

110-111 

111 

011-111 

t    T 

t    t 

Output 

Input 

D.   de  BRUIJN  DIAGRAMS 

The  de  Bruijn  diagram  compactly  contains  all  the  above 


n 


information.   The  diagram  contains  2   vertices  corresponding 


,n 


to  the  2   states  and  two  directed  edges  from  each  state  to 
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the  possible  pair  of  successor  vertices.   As  a  result  two 
arrows  exit  every  vertex  leading  to  conjugate  states,  and 
two  arrows  enter  every  vertex  arriving  from  companion  states 
For  an  example  consider  the  de  Bruijn  diagram  G   for  n  =  3 
in  Fig.  2.4.   An  Eulerian  path  through  G   is  defined  as  a 


path  which  visits  each  edge  exactly  once.   A  path  which 
visits  every  vertex  of  G   is  called  a  Hamiltonian  path.   A 
Hamiltonian  path  in  G   is  a  de  Bruijn  sequence  of  span  n. 
If  the  edges  of  G   are  labeled  with  an  n  +  1-tuple  defined 
by  the  labels  on  the  predecessor  and  successor  states,  then 
the  Eulerian  paths  in  G   correspond  to  Hamiltonian  paths- - 
and  de  Bruijn  sequences  in  the  graph  G  +, .   The  Hamiltonian 
path  in  Fig.  2.5  is  a  subgraph  of  G-  in  Fig.  2.4  with  exactly 
one  edge  emanating  from  each  vertex. 

E.   OPERATORS 

For  S  =  s,  s~  •••  s   ,  s   where  s.  e  {0,1}  for  l<i<n  the 
12      n-ln        l  —  — 

following  operators  are  defined. 

1 .  Identity  Operator 

The  identity  operator  (e)  operating  on  S  denoted 

eS  =  s,  s9  • • •  s   ,  s   =  S. 
12      n-ln 

2 .  Reverse  Operator 

The  reverse  operator  (r)  operating  on  S  is  denoted 

r  S,  where  r  S  =  s   s   ,  •••  s0  s, .   The  reverse  operator 

n  n- 1      2   1  r 

preserves  the  weight,  parity  and  length  of  S.   Note  also 
that  (r)2  =  e. 

3 .  Complement  Operator 

The  complement  (dual)  operator  on  S  is  denoted  d  S 

or  S,  where  S  =  s",  s -  •  •  •  S      ,  ?   such  that  s.  ■  s-  @  1 

12  n-ln  11 

(mod   2)    for   l<i<n. 
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Fig.  2.5   A  Hamiltonian  Path  through  G, 
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Lemma  1:   If  W  (S)  =  k  for  S  =  s,  s~  •••  s   then 

'  1   Z      n 

W  (S)  =  n-k. 

Proof:   W  (S)  =  k  implies  k  nonzero  elements  £  S 

and  n-k  zero  elements  e  S. 
Complementing  n-k  zero  elements  produces 

n-k  nonzero  elements  e  S". 
: .  W  (S)  =  n-k 

Q.E.D. 

Lemma  2:   If  P  (S)  =  a  where  a  e  {0,1}  and  S  =  {s,  s~  •••  s  }, 

1  l  n  7 

then  P  (S)  is  a  iff  n  is  even, 
n 
Proof:   P  (S)  =   E   s.  (mod  2)  =  a 

i  =  l   x 

n  n 

P  (3")  =   E   s.  (mod  2)=   E   s.  0  1  =  a  @  n 
i=l  i=l  1 

therefore  P  (5)  =  P  (n-a) 

Q.E.D. 

2 

Also  note  (d)   =  e,  and  that  the  reverse  and  comple- 
ment operators  are  commutative,  i.e.  rd  =  dr. 
4 .   Reverse  Complement  Operator 

The  reverse  complement  operator  on  S  is  denoted 
r  S  =  s"  s   -  •••  s~2  s"i  •   Weight  and  parity  are  effected 
by  the  reverse  complement  operator  exactly  as  with  the 
complementation  operator. 

These  operators  apply  to  sequences  as  well,  since 
sequences  are  composed  of  n-tuples.   Table  2.6  shows  the 
effect  of  operators  on  de  Bruijn  sequences  of  span  5  and  4 
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where  r  S,  S,  and  r  S  are  all  de  Bruijn  sequences,  as  they 
always  are  if  S  is  a  de  Bruijn  sequence. 

TABLE  2.5 
de  BRUIJN  SEQUENCE  OPERATORS 

n  =  3  (Canonical  Form) 

S  =  00011101  =        00011101 

rS  =  10111000  =        00010111 

S  =  11100010  =        00010111 

rS  =  01000111  =        00011101 

n  =  4  (Canonical  Form) 

S  =  0000111101001011  =  0000111101001011 

rS  =  1101001011110000  =  0000110100101111 

S  =  1111000010110100  =  0000101101001111 

rS"  =  0010110100001111  =  0000111100101101 

If  as  in  the  case  for  n  =  3  in  Table  2.5  a  sequence 
S  is  equivalent  to  its  reverse  complement,  then  S  is  termed 
a  RC  sequence.   This  applies  to  both  de  Bruijn  and  non- 
de  Bruijn  sequences.   The  following  theorem  appears  in  a 
paper  by  Etzion  and  Lempel  [Ref.  6]. 

Theorem  1:   A  sequence  S  is  a  RC  sequence  if  I    (S)  is 
even  and  S  =  {X,  rX}  for  some  X  where  X  is  a  binary  string. 
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F.   LINEAR  RECURSION  FORMULA 

de  Bruijn  sequences  generated  by  a  nonlinear  function 

can  also  be  generated  by  a  linear  function  f  (x,  • • •  x  )  =  x 
&  j  ^  1      n^     o 

having  the  linear  recursion  formula 

n 

S   a.  x.  =  0  . 

th 
Here  each  a.  determines  a  tap  to  the  j  L   stage  of  the 

register.   The  linear  generator  polynomial  must  be 

(x  +  1)   for  some  positive  integer  k,  due  to  the  periods  of 

de  Bruijn  sequences  (2  )  and  results  from  theorems  given 

by  Golomb  [Ref.  7:  pp.  27-43]. 

1 .   Pascal  Triangle  (Mod  2) 

The  values  for  each  a-  are  summarized  in  the  Pascal 

triangle  (mod  2)  Table  2.6  for  various  values  of  n.   They 

are  the  coefficients  generated  by  the  binomial  expansion  of 

(x  +  l)k  mod  2. 


TABLE  2.6 
PASCAL  TRIANGLE  (MOD  2) 

k  =   0  1 

1  11 

2  10         1 

3  1111 

4  10        0        0        1 

5  110        0         11  (x+1)5 

ao      al      a2      a3      a4      a5 
26 


2 .   Pascal  Triangle  Properties 

(a)  The  a-  sequence  is  palindromic,  that  is 
aQ  =  an  =  1,  a1  =  an_1,  a2  =  aR_2,  etc. 

(b)  The  number  of  nonzero  coefficients  is  even, 
n 

that  is   Z   a.  =  0  (mod  2) . 
j=0   J 
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III.   COMPLEXITY  OF  de  BRUIJN  SEQUENCES 

A.  DEFINITION  OF  COMPLEXITY 

The  term  complexity  indicates  a  general  measure  of  the 
predictability  of  a  sequence.   Various  definitions  of  com- 
plexity are  in  use,  but  this  paper  will  use  the  one  given 
by  Chan  and  Games  [Ref.  1].   Their  definition  of  complexity 
is  the  length  of  the  shortest  LFSR  required  to  generate  a 
sequence.   They  apply  their  results  primarily  to  de  Bruijn 
sequences.   This  definition  is  also  the  one  used  by  Herlestom 
[Ref.  8]  to  look  at  two  different  shift  register  cycle  gen- 
erators.  For  a  de  Bruijn  sequence  S,  let  C(S)  denote  the 
complexity  of  S. 

Previous  results  by  Chan  and  Games  [Ref.  1]  establish 
lower  and  upper  bounds  of  2    +n  and  2    respectively  for 
the  complexity  of  a  de  Bruijn  sequence.   Thus,  the  length 
of  a  LFSR  is  almost  as  long  the  sequence  as  opposed  to  a 
n-stage  NFSR  which  generates  the  sequence.   Table  3.1  serves 
as  a  ready  reference  for  the  upper  and  lower  limits  of 
complexity  for  3<n< 8 .   The  upper  bound  is  known  to  be  attained 
for  all  n  and  the  lower  bound  for  all  n<_6 . 

B.  COMPLEXITY  ALGORITHM 

A  fast  method  for  determining  the  complexity  of  a  de 
Bruijn  sequence  was  developed  by  Games  and  Chan  [Ref.  9] 
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TABLE  3.1 
LIMITS  OF  COMPLEXITY 

Span  Lower  Bound  Upper  Bound 

(n)  (2n'1+n)  =  (2n-l) 

3  7  7 

4  12  15 

5  21  31 

6  38  63 

7  71  127 

8  136        •  255 

and  is  represented  by  the  flow  chart  in  Fig.  3.1.   Let  S  be 

a  de  Bruijn  sequence  then  £(S)  =  2n.   Further  let  S  =  A:B  =  D 

where  A  =  {a,  a0  •  •  •  a   , }  and  B  =  {bn  b-,  •••  b    -, )  then 
1  l  ^n-1  1  L  ? 

D  =  (a,  a0  • • •  a   -,  b,  b->  •••  b    , }.   To  serve  as  a  guide 
12      ^n-112      ^n-1  ° 

an  example  for  n  =  5  is  now  presented.   Let  the  de  Bruijn 
sequence  S  =  {11111000001000110101100101001110}  having  length 
I    =  32. 

A   1111100000100011     I    =  16 

B    0101100101001110     C  =   0 

A+B   1010000101101101     C  =  16 

A   10100001  I    =   8 

B    01101101  C  =  16 

A+B   11001100  C  =  24 
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A 

1100 

B 

1100 

A+B 

0000 

A 

11 

B 

00 

I    =   4       Since  A©B  =  0 

C  =  24 

C  =  24       C  is  unchanged 


I    =   2 
C  =  24 
A+B   11  C  =  26 


A 

1 

B 

1 

A+B 

0 

i    =      1       Add  1  since  1    =    1 
C  =  26       5  A  =  1 
C  =  26       +  1  =  27 


C.  COMPLEXITY  DISTRIBUTION 

The  nature  of  the  distribution  of  complexity  of  de  Bruijn 
sequences  is  the  primary  interest  of  this  paper.   Let 
a  (c,n)  denote  the  number  of  de  Bruijn  sequences  having  com- 
plexity c  and  span  n.   The  complexity  distribution  of  de 
Bruijn  sequences  of  span  n  for  3<_n<_6  are  listed  in  Table 

3.2.   Unfortunately  the  large  number  of  de  Bruijn  sequences 

57 
for  n>7   (2   )  does  not  allow  an  exhaustive  examination  of 

their  complexities. 

Examination  of  this  data  led  Chan  and  Games  to  conjecture 

that  for  n>_3  a  (c,n)  is  congruent  to  0  modulo  4. 

D.  EFFECTS  OF  OPERATORS  ON  COMPLEXITY 

Chan  et  al .  [Ref.  1]  proved  that  a(c,n)  for  n>_3  was  con- 
gruent to  0  modulo  2.   By  showing  that  for  each  de  Bruijn 
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D    =    S 

t    =    2n 
c    =    0 

D    =   A:B 

I    =  1/2 
B  =  A@B 


es 


c  =  c  +  1 


Fig.  3.1    Flowchart  for  Complexity  Algorithm 
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TABLE  3.2 
COMPLEXITY  DISTRIBUTION 


n  =  4 

n 

=  5 

c 

a(c,4) 

c 

a(c,5) 

c 

a(c,5) 

12 

4 

21 

9 

27 

64 

13 

0 

22 

0 

28 

180 

14 

4 

23 

12 

29 

224 

15 

8 

24 

20 

30 

448 

25 

32 

31 

1024 

26 

36 

n  =  6 
a (c ,6)         c      a (c ,6) 


38 

448 

51 

8704 

39 

o  - 

52 

18096 

40 

32 

53 

34224 

41 

96 

54 

67700 

42 

160 

55 

126592 

43 

80 

56 

259320 

44 

432 

57 

519752 

45 

288 

58 

1041252 

46 

896 

59 

2090716 

47 

1168 

60 

4162352 

48 

2772 

61 

8342176 

49 

2352 

62 

16692832 

50 

5224 

63 

33731200 
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sequence  having  complexity  c,  there  exists  another  different 

de  Bruijn  sequence  S  also  having  complexity  c.   Consequently 

they  showed  that  the  complement  operator  does  not  change  the 

complexity  of  a  de  Bruijn  sequence.   Using  a  different 

approach,  the  next  theorem  will  show  that  the  complexity  of 

a  de  Bruijn  sequence  is  also  the  complexity  of  the  de  Bruijn 

sequences  defined  by  the  reverse  operator  and  by  the  reverse 

complement  operator  as  well. 

Lemma  3 .   For  S  =  {s,  s2  •••  s   _ , }  ,  a  de  Bruijn  sequence, 

c  2 

C(S)  =  c  iff  E    a-  s.  =  0  for  everv  s.  where  each  a-  comes 
-nil  i  i 

th  1  = 
from  the  c   row  of  the  Pascal  triangle  (mod  2)  and  c  is  the 

smallest  integer  such  that  the  above  recursion  holds. 

Proof.   S  is  de  Bruijn  with  C(S)  =  c  iff  s.  =  f(s-  ,••• 
J        *■  '  l     v  l  +  l 

S-.  )  for  every  s.  e  S,  which  is  equivalent  to  s  -  = 
l+c  J       i  ^  1 

c 

E   a.  s.    from  results  by  Golomb  [Ref.  5:  pp.  27-43] 

j=l  ^   1+J  c 

mentioned  earlier.   In  turn  s.  =  E   a-  s-  •  =  0  from  proper- 

1   j=0  3      1   J 

ties  of  the  Pascal  triangle  given  in  Section  II.     Q.E.D. 

This  lemma  is  used  in  the  following  theorem  to  show  that 
complexity  is  preserved  by  various  operators  on  de  Bruijn 
sequences . 

Theorem  2.   If  S  is  a  de  Bruijn  sequence  having  complexity 
C,  then  C(S)  =  C(5)  =  C(rS)  =  C(rS")  =  c. 

Proof.   Part  I  (reverse  operator) 

Let  s.  s.+,  •••  s-    be  an  arbitrary  (c+1)  long  segment 
of  S  which  satisfies 
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c 
(1)   S   a-  s.  .    n 

j=0   J   i+J  =  o 

by  Lemma  3.   Since  the  a-  sequence  is  palindromic  as  shown 
earlier,  s-    •••  s-  ,  s.  a  (c+1)  long  segment  of  rS  satisfies 

l+C        1+1   1 

equation  (1).   Thus  C  (rS)  =  c. 

Part  II  (complement  operator)  substituting  s".  •••  s".+ 
the  c  +  1  long  segment  of  S",  into  equation  (1)  ,  the  results  are 

ao  ?i  +  al  ?i  +  l  +  ••'  +  Vl  ?i+c-l  +  ac  ?i  +  c  = 

aQ  (si+  1)  ♦  ax  (si+1  +  1)  *  ...  +  ac_1  (si+c-1  *1)  +  ac 

c  c 

(s..„  +1)  =   Z   a-  s.  .  +  £     a-  (mod  2).   Since  there  are 

J  =0    J       J     J  =0    J 

C 

an  even  number  of  nonzero  a-s,   E   a.  =  0  and  equation  (1) 

J    j=0   3 
is  satisfied.   Therefore,  by  Lemma  3  C(£)  =  c. 

Part  III  (reverse  complement  operator) ,  by  Parts  I  and 

II,  it  follows  immediately  that  C  (rS")  =  c.   Thus  C  (S)  = 

C  (rS)  -  C  (?)  =  C  (rS) 

Q.E.D. 
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IV.   TRUTH  TABLES 

Linear  and  nonlinear  feedback  functions  were  used  in 

conjunction  with  n-stage  shift  registers  to  generate  de  Bruijn 

sequences  in  Section  II.   The  input  xfi  was  generated  as  some 

function  f  of  the  current  state  x,  x0  •••  x  ,  that  is 

11  n ' 

(1)  xQ  =  f  Cx1  x2  •••  xn) 

or  equivalently  writing  f  as  in  (2)  when  specializing  for 
cycles  only  functions. 

(2)  xQ  =  xn  ♦  g  (Xl  •••  xn.1). 

Extensive  use  is  made  of  equation  (2)  since  results  for 
cycles  only  are  desired. 

A.   FULL  TRUTH  TABLE 

The  full  truth  table  is  a  useful  way  of  listing  the 
functional  value  x~  from  equation  (1)  for  each  possible 
state  of  the  shift  register.   The  functional  values  in  Table 
4.1  is  an  example  of  a  full  truth  table  for  a  de  Bruijn 
sequence  (00011101) . 

Examining  the  full  truth  table  it  becomes  apparent  that 
the  information  in  the  lower  half  (below  the  dashed  line) 
is  redundant,  since  the  values  of  f  (x,  x~  x,)  for  x,  =  1 
in  the  lower  half  are  merely  the  complements  of  the  values 
of  f  (x,  Xy   xJ  for  x,  =  0  in  the  upper  half  for  the  identical 
values  of  x,  x2 .   This  suggests  that  the  upper  half  (or  the 
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TABLE  4.1 
FULL  TRUTH  TABLE  OF  (00011101) 

Current  State       Functional  Value 
Cx5  x2  xl} xQ  =  f  (x1  x2  x5) 

0   0   0  1 

0   0   1  1 

0   10  0 

_0_  1  _1_  _1_ 

10   0  0 

10   1  0 

110  1 

111  0 

lower  half)  of  the  truth  table  can  display  all  the  informa- 
tion for  a  function  generating  only  cycles. 

B.   HALF  TRUTH  TABLE 

Focusing  on  equation  (2) ,  since  x   is  always  1  in  the 
lower  half  this  can  be  viewed  as  complementing  the  functional 
values  for  g  (x,  •••  x   ,)  in  the  upper  half.   In  this  manner 
the  half  truth  table  in  Table  4.2  is  constructed  for  the  same 
de  Bruijn  sequence  (00011101)  to  aid  comparison  with  Table 
4.1.   From  this  point  on  the  term  truth  table  will  always 
mean  the  upper  half  truth  table  with  x   =0,  unless  noted 
otherwise . 
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TABLE  4.2 
HALF  TRUTH  TABLE  OF  (00011101) 


x2 

xl 

cr 
a 

(Xj    x2) 

0 

0 

1 

0 

1 

1 

1 

0 

0 

1 

1 

1 

The  string  of  functional  values  in  the  truth  table  is 


defined  to  be  the  generator   G  =  {gn  g,  •••  g 


,n-l 


}  where 


211  *-l 

the  subscripts  of  g  are  the  decimal  equivalent  of 

x   ,  •••  x,  .   Both  the  weight  and  parity  of  generator  G 
will  be  of  interest  as  the  investigation  of  this  paper  con- 
tinues.  Whenever  the  parity  or  weight  of  the  truth  table 
is  mentioned  in  this  paper,  it  is  referring  to  the  generator 
of  that  truth  table. 

C.   ZERO  TRUTH  TABLE 

The  zero  truth  table  is  defined  as  the  truth  table  in 
which  all  functional  values  are  zero.   The  zero  truth  table 
for  n  =  3  contained  in  Table  4.3  is  basically  the  output  of 
a  circulating  register  which  generates  pure  cycles. 

Pure  cycles  are  defined  as  cyclic  permutations  of  the 
original  state.   The  zero  truth  table  in  Table  4.3  generates 
the  four  pure  cycles  shown  in  Fig.  4.1.   The  number  of  pure 
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TABLE  4.3 
ZERO  TRUTH  TABLE  FOR  n  =  3 

x2  x1  g  (x1   x2) 

0   0  0 

0   1  0 

10  0 

11  0 

cycles  Z   has  been  shown  [Ref.  7:  p.  120]  to  be 

Z   =  1   2   <D  (d)  2   a. 
n  d|  n 

Eulers  totient  function  $  (d)  is  the  number  of  fractions 
of  the  form  a/d  where  lfa<d  and  a/d  is  in  lowest  terms.   The 
summation  is  over  all  positive  "integers  d  which  divide  n, 
denoted  (din).   It  is  further  shown  [Ref.  5]  that  Z  -1  is 
the  minimum  possible  weight  truth  table  that  can  generate  a 
de  Bruijn  sequence. 

In  the  next  section,  the  proof  of  theorem  6  gives  an 

example  of  how  pure  cycles  can  be  joined  in  a  way  to  create 

a  de  Bruijn  sequence  that  has  a  minimum  weight  truth  table. 

Table  4.5  gives  a  listing  of  Z   and  Z  -1  for  l<n<7. 
to  to     n      n        —  — 

D.   ONES  TRUTH  TABLE 

The  ones  truth  table  assigns  a  functional  value  of  1  to 
every  entry  of  g  (x,  •••  x  , ) .   The  ones  truth  table  depicts 
the  output  of  a  complementing  shift  register.   The  input  is 
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Fig.    4.1      Pure   Cycles    (n=3) 
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always  the  complement  of  the  output,  i.e.  f  (x,  •••  x  )  = 

1  +  x  .   The  ones  truth  table  for  n  =  3  is  shown  in  Table 
n 

4.4.   The  complementing  cycles  from  Table  4.4  are  depicted 
in  Fig.  4.2. 

TABLE  4.4 
ONES  TRUTH  TABLE  (n  =  3) 


x2    xl 

a 

(xx   x2) 

0      0 

1 

0      1 

1 

1      0 

1 

1      1 

1 

The  number  of  complementing  cycles  Z   .[Ref.  7]  is 


n 


=  7   -  1 


$  (2d)  2 


ia" 


2n  2d  n 


Since  summation  for  Z    is  over  all  2d  (even)  numbers  which 

n 


divide  n,  for  n  odd  Z 


n 


■  V 


Accordingly  Fredricksen  [Ref.  5]  shows  truth  tables  of  max 

weight  that  can  generate  a  de  Bruijn  sequence  is 

211-1  -  Z  *  +  1. 
n 


The  values  of  Z   and  2n_1 

n 


Z   +1  are  given  in  Table  4.5 
n  ° 


for  l£n<7.  Also  included  are  the  number  of  de  Bruijn  se 
quences  having  truth  tables  of  maximum  or  minimum  weight 
from  a  listing  by  Fredricksen  [Ref.  5]. 
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TABLE  4.5 
PURE  CYCLE  COMPILATION 


#  de  Bruijn 

#  de  Bruijn 

Sequences 

Sequences 

n 

Zn 

Z  -1 
n 

with  Min 
Wt  Truth  Table 

7  * 

n 

2n"1-z  *+i 
n 

with  Max 
Wt  Truth  Table 

1 

2 

1 

1 

1 

1 

1 

2 

3 

2 

1 

1 

2 

1 

3 

4 

3 

2 

2 

3 

2 

4 

6 

5 

12 

2 

7 

3 

5 

8 

7 

26-32 

4 

13 

26 

6 

14 

13 

,14  ,4  .2 
2   •  j  •  5 

6 

27 

214' 

7 

20 

19 

228«35-53«13 

10 

55 
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Fig.  4.2   Complementing  Cycles  (n=3) 

E.   EFFECTS  OF  CHANGING  THE  TRUTH  TABLE 

If  the  ones  truth  table  in  Table  4.4  is  changed  so  that 
g  (0,1)  =  0  vice  1,  the  truth  table  of  Table  4.6  is  produced 
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TABLE  4.6 
TRUTH  TABLE  (TABLE  4.4  MODIFIED) 


x2    xl 

g 

(x1    x2) 

0      0 

1 

0      1 

0 

1      0 

1 

1      1 

1 

This  change  in  the  truth  table  effects  the  cycles  of  Fig. 
4.2  by  joining  them  as  depicted  in  Fig.  4.3.   Indeed  the 
change  produces  a  Hamiltonian  path  for  the  de  Bruijn  sequence 
(00010111)  . 

Suppose  the  truth  table  of  Table  4.6  is  now  changed  so 
that  g  (00)  =  0  vice  1.   This  produces  the  truth  table  in 
Table  4.7,  and  its  cyclic  composition  in  Fig.  4.4. 

TABLE  4.7 
TRUTH  TABLE  (TABLE  4.6  MODIFIED) 


x2    xl 

o 
o 

Cxx   x2) 

0      0 

0 

0      1 

0 

1      0 

1 

1      1 

1 
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Fig.  4.3   Cyclic  Composition  of  Table  4.4 
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Fig.  4.4   Cyclic  Composition  of  Table  4.7 
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This  change  in  the  truth  table  causes  the  cycle  to  split 
into  two  cycles.  Suppose  (x  •••  x,)  is  on  element  on  upper 
of  two  cycles  as  shown  in  Fig.  4.5  and  that  its  successor  is 


n-1 


Fig.  4.5   Two  Cycle  Graph. 


x,  1  as  f  (x. 


x  )  =  1.   If  the  truth  table  is 
n 


now  changed  so  that  f  (x,  •••  x  )  =0,  then  x   ,  •••  x,  0 
°  k  1      nJ  n-1      1 

is  the  new  successor  for  x   • • •  xn  .   Retaining  the  cycles 

n      l  ° 

only  condition  f  (x,  •••  x   ,  x  )must  now  be  1 .   So  the 
'  1      n  - 1   n 

successor  of  x  x   .  •••  x,  is  now  x   ,  •••  x-,  1  .   One  of 

n  n-1      1         n-1      1 

two  things  can  happen  how  depending  on  whether  x   ,  •••  x,  0 
is  on  the  upper  or  lower  cycle. 
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First  if  (x   ,  •••  x,0)  is  on  the  upper  cycle,  then  that 
cycle  is  split  as  shown  in  Fig.  4.6  forming  three  cycles. 


Fig.  4.6   Cycle  Splitting. 

Otherwise  if  (x   ,  •••  x,0)  is  on  the  lower  cycle,  then 

n- 1      1  J  ' 

those  two  cycles  are  joined  as  shown  in  Fig.  4.7. 

Clearly  the  parity  of  the  cycles  changes  for  each  change 
in  the  truth  table.   The  following  theorem  from  Golomb 
[Ref.  7:  p.  122]  shows  that  the  parity  of  the  number  of  cycles 
and  the  parity  of  the  truth  table  generator  are  equal. 

Theorem  3.   For  n>2,  the  parity  of  the  truth  table  is 
equal  to  the  parity  of  the  number  of  cycles  of  the  truth 
table. 
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Since  a  de  Bruijn  sequence  is  one  cycle  by  definition, 
the  following  corollary  of  theorem  3  is  given. 

Corollary.   If  a  sequence  is  de  Bruijn,  then  the  parity 
of  the  TT  of  that  sequence  is  odd. 


Fig.  4.7   Cycle  Joining. 

F.   SPECIAL  NUMBERS 

Analyzing  the  effects  of  de  Bruijn  sequence  operators 
on  their  truth  tables,  certain  binary  numbers  are  found  to 
have  special  properties  which  can  be  exploited  in  later 
analysis . 


48 


1.  Type  RC  Numbers 

If  X  =  x,  •••  x   is  a  binary  number,  such  that 
In 

x   •••x,   =  x-,  •  •  •  x   then  X  is  called  a  RC  number, 
n      1     1      n  n 

J 

Lemma  4 .   For  n  even  there  are  exactly  2"  type  RC  numbers 

Proof.   For  n  even,  if  X  is  an  RC  number,  then 

x,  • • •  x   =  x   • • •  x,  by  definition.   Therefore, 

x-i  =  x  ,  x~  =  x   ,  ,  etc.   Thus  x.,«*»x   can  be  written  as 
1    n'   2    n- 1 '  1    n 

x,x~»««x  x  •••x%x,  .   Since  ^  positions  can  be  filled  in 

iz    nn    z  l  z  r 

2   I 

either  of  two  ways,  the  total  number  of  RC  numbers  for  n 
n      J     ' 

-2 
even  is  2  . 

Q.E.D. 

There  are  no  type  RC  numbers  for  n  odd,  since  this 

would  require  t>hat  x   -  =  x  .-,  which  is  impossible. 
M  n+1    n+1  v 

2 .  Type  R  Numbers 

If  X  =  x^'-'X   is  a  binary  number  such  that 
1    n 


=  x 


•••x,,  then  it  is  called  a  type  R  number 
n    1 '  /r   n 


2 
Lemma  5 .   For  n  even  there  areexactly  2   type  R  numbers, 

— 2 

and  for  n  odd  there  are  exactly  2     type  R  numbers. 

Proof.   The  proof  for  n  even  parallels  that  of  Lemma  4 

and  will  not  be  repeated.   For  n  odd,  since  x  _,,  =  x  ^, 

^  '        n+1     n+1 

~^~         ~T~ 
there  are  n+1  positions  to  be  filled  in  either~pf . two 

ways.   Therefore,  for  n  odd  there  are  exactly  2  w   type  R 

numbers . 

Q.E.D. 
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A  useful  way  of  organizing  special  numbers  and  other 
truth  table  values  to  assist  in  evaluating  effects  on  the 
generators  for  various  operators  is  shown  in  Table  4.8  for 
n  =  5  and  n  =  6.   The  numbers  are  the  decimal  equivalent  of 
the  binary  (n-1) -tuple,  xn-i*,,x]>  from  the  feedback  function 

x0  =  xn  +  §  (xl  ••'  xn-l) 
which  also  corresponds  to  the  truth  table  state. 

For  example  for  n  =  5,  if  f  (2)  =  1  for  some  sequence 

S  then:   f  (11)  =  1  for  r  S,  f  (4)  =  1  for  r  S,  and  f  (13)  =  1 

for  S.   In  general,  for  some  operator  8  and  some  number  a, 

if  f  (a)  =  1  for  some  sequence  S,  then  f  (6a)  =  1  for  8  S. 

Special  numbers  have  properties  in  addition  to  this.   If  a 

is  a  type  RC  number  and  f  (a)  =  1  for  S,  then  f  (a)=  f  (r  a)  = 

1  for  r  S.   Similarly  if  y  is  a  type  R  number  and  f  (y)  =  1 

for  S,  then  f  (y)  =  f  (ry)  =  1  for  r  S.   Thus,  if  we  want  to 

shown  S  is  a  sequence  such  that  S  =  r  S,  then  if  f  (a)  =  1 

in  S,  then  f  (r  a)  =  1  in  S  also. 

G.   GENERATOR  ANALYSIS 

Analysis  of  the  generators  for  the  de  Bruijn  sequences 
S,  r  S,  S,  §  r  S  in  Table  4.9  appears  to  indicate  the  gener- 
ator G  of  sequence  S  denoted  G  (S~)  is  the  reverse  of  the 
generator  G  (S) .   Notice  also  that  G  (r  S)  =  r  [G  (r  S ) ] 
for  S.   Although  Table  4.9  is  an  example  of  a  single  sequence 
S,  it  is  not  an  isolated  incident;  in  general  the  result 
G  (S)  =  r  [G  (S~)  ]  holds  and  is  proven  in  the  following 
theorem. 
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TABLE  4.8 
TRUTH  TABLE  ANALYSIS  CHART 


s 

S 

rS 

rS 

s 

S 

rS 

rS 

rS 

rS 

S 

S 

rS 

rS 

S 

S 

0 

15 

0 

15 

1 

14 

8 

7 

2 

13 

4 

11 

3 

12 

12 

3 

5 

10 

10 

5 

6 

9 

6 

9 

0 

31 

0 

31 

1 

30 

16 

15 

2 

29 

8 

23 

3 

28 

24 

7 

4 

27 

4 

27 

5 

26 

20 

11 

6 

25 

12 

19 

9 

22 

18 

13 

10 

21 

10 

21 

14 

17 

14 

17 

Truth  table  states  (n=5) 
Type  R  numbers 
(0,6,9,15) 

Type  RC  numbers 
(3,5,10,12) 

*•  Truth  table  states  (n=6) 
Type  R  numbers 
(0,4,10,31) 
Type  RC  numbers 
(None) 
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TABLE  4.9 

GENERATOR  ANALYSIS  FOR  de  BRUIJN  SEQUENCES 

(r  S  f    S) 

S  11111001000111011010011000001010 

S  11111010100000110111000100101100 

rS  11111010100000110010110111000100 

rS  11111001101001000111011000001010 


X 

n- 

•  •  •    x  _ 

1,         ,xl 

G    (S) 

G    (S) 

G    (rS) 

G    CrS") 

0 

1 

1 

1 

1 

1 

0 

1 

1 

0 

2 

1 

0 

0 

1 

□ 

S 

1 

0 

0 

1 

Type   RC 

4 

0 

1 

1 

0 

lh 

0 

1 

1 

0 

O 

© 

0 

1 

0 

1 

Type   R 

7 

0 

1 

1 

0 

/other  than\ 

|        8 

1 

0 

0 

1 

\0  §  15    S 

© 

1 

0 

1 

0 

1,10,1 

1 

0 

0 

1 

11 

1 

0 

0 

1 

|l2| 

0 

1 

1 

0 

13 

0 

1 

1 

0 

14 

1 

0 

0 

1 

15 

1 

1 

1 

1 
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Theorem  4:   If  S  is  a  de  Bruijn  sequence  with  generator 
G  (S) ,  then  S  is  a  de  Bruijn  sequence  with  generator 
r  [G  (S)]. 

Proof:   Let  G  (S)  =  {g,,  g2,  •••»  g  n_i^   then  the  full 
truth  table  for  S  is  given  in  Table  4.10. 


TABLE  4.10 
FULL  TRUTH  TABLE  FOR  S 


n   n-1   n-2 


f  (xr, 


^ 


§1 
§2 


«2n-l 


§2 


a 

*2n-l 


If  x  ,  x   ,  , 
n'   n-1 ' 


,  xx,  f  (x1 


x  )  is  a  (n+1)  long  string 


in  S,    then  x    ,    x      ,,    •••,    x,  ,    T    Tx~I    xT    TTl    x   )      is   a    (n+1) 
n*      n-1  1  ^    1      2  nJ  v        ' 

long  string  in  S". 

Therefore  we  can  transform  each  line  of  the  full  truth 
table  for  S  by  taking  its  complement  to  get  Table  4.11. 
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TABLE  4.11 
FULL  TRUTH  TABLE  FOR  S  (NON  STANDARD) 


X 

n 

X 

n- 

•1 

X 

n- 

•2 

•   •   • 

xl 

f    Cx-p    x2,    •• 

•    ,   x   ) 

1 

1 

•   •    • 

1 

*1 

1 

1 

1 

•    •    • 

0 

h 

• 

•   • 

• 

• 

•   • 

• 

• 

•   • 

• 

1 

0 

0 

•    •   • 

0 

V"1 

0 

1 

1 

•   •   • 

1 

*i 

0 

1 

1 

•   •   • 

0 

*2 

• 

•    • 

• 

• 

•   • 

• 

• 

•   • 

• 

0 

0 

0 

•   •   • 

0 

a 

«2n-l 

Rearranging  the  rows  from  top  to  bottom,  the  TT  for  S  is 
placed  in  standard  form  as  depicted  in  the  Table  4.12. 
Further  examination  of  the  generators  in  Table  4.9  shows 
the  structure  of  G  (S)  to  be  remarkably  similar  to  the 
structure  of  G  (r  S) ,  with  the  exception  that  the  functional 
values  of  type  R  numbers  (6  5  9)  are  transposed.   If,  however, 
these  values  should  be  the  same  then  clearly  G  (5~)  =  G  (r  S) 
or  S  =  r  S.   Indeed,  this  occurs  twice  for  n  =  3  and  64  times 
for  n  =  5,  and  will  be  shown  to  happen  for  all  odd  n> 3 . 
This  cannot  happen  for  n  even,  because  there  are  no  type 
RC  numbers  in  the  truth  table. 


54 


TABLE  4.12 
FULL  TRUTH  TABLE  FOR  S 


X 

n 

V 

•1 

X 

n- 

-2 

•   •   • 

xl 

f    (xx,    x2,    " 

•     ,    x    ) 

0 

0 

0 

•   •   • 

0 

g7n-l 

• 

• 

» 

• 

• 

• 

• 

• 

• 

• 

• 

• 

0 

1 

1 

•    •   • 

0 

§2 

0 

1 

1 

•   •   • 

1 

*1 

1 

• 
• 

0 

0 

•  • 

•  • 

•   •   • 

0 

a 

%n-l 

• 

• 

•   • 

• 

1 

1 

1 

•    •   • 

0 

§2 

1 

1 

1 

•   •    • 

1 

Si 

Clearly  G  (S)  is  r  [G  (S) ] 


Q.E.D 
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V.   RESULTS 

Examining  Table  3.2  it  can  be  quickly  determined  that  the 
number  of  de  Bruijn  sequences  of  given  complexity  for  3<n<_6 
is  0  (mod  4).   Games  and  Chan  [Ref.  1]  made  this  conjecture, 
but  were  able  to  show  only  that  the  distribution  a  (c,  n) 
was  0  mod  2.   Their  result  followed  by  showing  that  for  a 
given  de  Bruijn  sequence  S  of  complexity  c,  there  existed 
a  unique  de  Bruijn  sequence  S  also  having  complexity  c. 
In  Section  III  it  was  shown  that  reverse  and  reverse  dual 
operators  preserve  the  complexity  of  de  Bruijn  sequences  as 
well. 

If  each  of  these  new  operators  produced  unique  de  Bruijn 
sequences,  then  it  would  be  easy  to  show  that  the  numbers 
a  (c,  n)  e  0  (mod  4)  since  for  every  de  Bruijn  sequence  S 
there  would  also  exist  S,  r  S  and  r  S  all  de  Bruijn  and 
unique  having  complexity  c.   As  mentioned  in  Section  IV,  the 
sequences  r  £  and  r  S  are  not  always  distinct  from  the 
sequences  S  §  J  respectively;  however,  for  n  even  the  next 
theorem  shows  that  all  four  of  these  sequences  are  distinct. 
In  addition  it  will  be  shown  that  RC  sequences  occur  for 
all  n  odd,  and  a  lower  bound  for  the  number  of  RC  sequences 
will  be  established  for  all  n  and  investigated  for  n  =  7. 
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A.   a  (c,  n)  =    0  MODULO  4  FOR  n  EVEN 

Etzion  and  Lempel  [Ref.  6]  give  a  proof  that  a  (c,  n) 
is  congruent  to  0  mod  4  for  c  n  even.  In  the  proof  for  n 
even  given  here  a  different  approach  is  taken. 

Theorem  5.   For  even  n>_ 3 ,  the  numbers  a  (c,  n)  e  0  (mod  4) 
for  a  de  Bruijn  sequence  S  of  span  n  and  complexity  c. 

Proof.   From  theorem  1,  it  is  known  that  for  a  de  Bruijn 
sequence  S  of  complexity  c,  that  S,  r  S  §  r  S  also  have 
complexity  c.   Games  and  Chan  showed  the  uniqueness  of  5", 
so  all  that  remains  is  to  show  that  r  S  is  distinct  from 
S,  or  equivalently  that  r  S  is  distinct  from  S  for  n  even. 

Suppose  on  the  contrary  that  S  =  r  S".   Then  consider 

an  n-tuple  Oa   ,  • • •  a,  in  the  truth  table  such  that 
r     n-1      1 

f  (a,  •••  a   ,  0)  =  1.   Then  0a   ,  •••  a,  1  is  a  n+1  long 
v  1      n-1   '  n-1      1  & 

string  in  S,  and  o  a,  •••  a   ,  1  is  a  n+1  long  string  in  r  3". 

If  S  =  r  S 

then  0  a,  • • •  a   , 1  is  a  n+1  long  string  in  S, 
1      n- 1  o  o 

which  implies 

that  I  a.  •••  a  -i  0  is  a  n+1  long  string  in  S,  and 

therefore  f  (a,  •••  a   .  0)  =  1.   This  must  hold  for  every 
1      n- 1 

x,  •••  x„  ,  0  such  that  f  (x,  •••  xM  ,  0)  =  1.   Pairing 
1      n-1  1      n-1  ° 

these  vectors,  a,  •  •  •  a   ,  and  a   ,•••!",,  two  at  a  time 

1      n-1      n-1      1 

results  in  an  even  number  of  l's  in  the  truth  table.   So 

there  must  be  a  vector  which  pairs  with  itself  as  the  parity 

of  ones  is  odd.   But  then  x   ,  • • •  x,  =  x,  • • •  x   ,  is  a 

n-1      1    1      n-1 
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contradiction  since  x  ^,  f   x  .,  for  n  even.   Therefore  for 

n+1    n+1 

2       2 
n  even  S  7^  r  S  and  r  S  f   S,  so  r  S  and  r  S  are  distinct 

and  a  (c,  n)  =0  mod  4. 

Q.E.D. 

B.   RC  SEQUENCES  EXIST  FOR  ALL  ODD  n>3 

The  results  a(c,n)  =    0  (mod  4)  for  all  n  are  not  obtained 
as  hoped  for.   For  n  =  3  and  n  =  5  RC  sequences  are  known 
to  exist,  but  it  is  not  known  if  this  is  true  for  all  odd 
n^3  or  just  some.   The  following  theorem  shows  that  RC 
sequences  exist  for  all  odd  n>_3 .   Recall  that  RC  sequences 
are  those  sequences  S  for  which  S  =  r  S  or  equivalent  by 
S  =  r  S  since  r  (S)  =  r  (r  3")  yields  S  =  r  S. 

Lemma  6 .   The  sequence  formed  by  joining  all  pure 
cycles  of  n-tuples  is  de  Bruijn. 

Proof.   The  least  elements  of  each  pure  cycle  for  a 
given  n  are  easily  composed  and  arranged  in  order  by  weight 
according  to  the  algorithm  in  Appendix  A.   Since  the  pure 
cycles  contain  all  n-tuples  exactly  once,  all  that  remains 
is  to  show  how  each  pure  cycle  can  be  interconnected,  thus 
forming  a  de  Bruijn  sequence.   For  each  cycle  of  wt  a>_l 
there  exists  a  least  element  "2X+1"  (necessarily  odd,  if 
it  were  even  (2X)  then  "XM  is  on  the  same  cycle  and  smaller) 
whose  predecessor  "X"  has  wt  (a-1) .   Since  the  pure  cycles 
are  ordered  by  weight,  it  is  clear  that  the  least  element 
of  each  pure  cycle  of  wt  a>_l  can  be  joined  to  a  cycle  of 
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wt  (a-1)  by  changing  the  truth  table  for  its  predecessor, 
i.e.  g(X)  becomes  1.  In  this  way  all  the  pure  cycles  are 
now  interconnected,  and  the  sequence  is  de  Bruijn.. 

Q.E.D. 

Theorem  6 .   For  all  odd  n>_3 ,  there  exist  RC  de  Bruijn 
sequences . 

Proof.   An  RC  de  Bruijn  sequence  can  always  be  composed 

in  the  following  way  for  n>_3 .   First  construct  a  light  cycle 

(LC)  for  those  pure  cycles  of  wt  a<_— j—        identical  to  the 

method  used  in  Lemma  6.   Note  that  when  the  cycles  of 

wt  —j—       are  joined  to  cycles  of  smaller  weight  the  truth 

table  is  changed  at  positions  whose  weight  is  less  than 

— y—   .    For  every  X  such  that  f  (X)  =  1  in  LC,  induce 

_  n+1 

f  (r  X)  =  1  on  the  pure  cycles  of  wt  (n-a)  >_  2   ,  forming 

a  heavy  cycle  (HC) .   Note  that  HC  is  the  reverse  complement 

of  LC.   All  of  the  positions  changed  so  that  g  (X)  =  1  on 

HC  are  of  weights  bigger  than  or  equal  to  — j—   • 

Join  LC  and  HC  by  using  a  type  RC  number,  y,  which  exist 

by  Lemma  4,  since  for  n  odd  the  truth  table  states  have 

even  (n-1)  length.   The  RC  number  y  has  weight  ^-j—   •   Letting 

f  (Y)  =  1  then  joins  LC  and  HC  forming  one  sequence.   By 

Theorem  1  the  sequence  is  RC,  and  by  Lemma  6  the  sequence 

is  de  Bruijn. 

Q.E.D. 
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C.   LOWER  BOUND  ON  NUMBER  OF  RC  SEQUENCES 

1 .  Theorem 

A  weak  lower  bound  on  the  number  of  RC  de  Bruijn 
sequences  is  established. 

Theorem  7:   A  lower  bound  on  the  number  of  de  Bruijn 
RC  sequences  for  odd  n>3  is  A  •  2  "   ,  where  A  is  the 
number  of  possible  different  interconnections  for  pure 
cycles  of  weighty— ^—  as  determined  by  the  Best  method 
[Ref.  5]. 

Proof:   Restricting  attention  to  those  cycles  of 
wt<_— j—      (since  the  other  connections  are  induced)  ,  the 
Best  method  gives  the  number  A  of  possible  different  inter- 
connections for, the  pure  cycles  of  weight<^y—  .Since 
there  exist  2     type  RC  numbers,  this  gives  2     ways 

of  joining  A  different  pairs  of  cycles  for  a  total  of 

2 
A  •  2     possibilities. 

Q.E.D. 

2 .  Example  for  n  =  7 

The  pure  cycles  of  wt<_3  and  the  adjacency  matrix 
of  A  are  given  in  Tables  5.3  and  5.4  respectively.   The 
value  of  the  determinant  of  the  adjacency  matrix  is  38,880 
for  any  diagonal  element.   This  gives  a  lower  bound  of 
38,880  X  8  =  311,040  for  the  number  of  de  Bruijn  RC  sequences 
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TABLE  5.1 
PURE  CYCLES  FOR  N=7  $  wt  <  3 


0 

1 

2 

4 

8 

16 

32 

64 

3 

6 

12 

24 

48 

96 

65 

5 

10 

20 

40 

80 

33 

66 

7 

14 

28 

56 

112 

97 

67 

9 

18 

36 

72 

17 

34 

68 

11 

22 

44 

88 

49 

98 

69 

13 

26 

52 

104 

81 

35 

70 

19 

38 

76 

25 

50 

100 

73 

21 

42 

84 

41 

82 

37 

74 
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TABLE 

5 

.2 

ADJACENCY 

MATRIX  FOR 

n=7 

1 

-1 

0 

0 

0 

0 

0 

0' 

0 

0 

1 

7 

-2 

-2 

-2 

0 

0 

0 

0 

0 

0 

-2 

7 

0 

0 

-2 

-1 

-1 

-1 

0 

0 

-2 

0 

7 

0 

-1 

-1 

-1 

0 

-2 

0 

-2 

0 

0 

7 

0 

-1 

-1 

-2 

-1 

0 

0 

-2 

-1 

0 

3 

0 

0 

0 

0 

0 

0 

-1 

-1 

-1 

0 

3 

0 

0 

0 

0 

0 

-1 

-1 

-1 

0 

0 

3 

0 

0 

0 

0 

-1 

0 

-2 

0 

0 

0 

3 

0 

0 

0 

0 

-2 

-1 

0 

0 

0 

0 

3 
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VI.   CONJECTURES 

Although  a  complete  proof  for  a  (c,  n)  =    0    (mod  4) 
has  yet  to  be  given,  a  technique  for  grouping  RC  sequences 
in  a  0  (mod  4)  fashion  will  be  demonstrated  for  n  =  5. 
However,  the  technique  fails  to  preserve  complexity  in 
every  case.   However,  for  n  =  5  the  complexitites  produced 
were  congruent  to  0  (mod  2)  and  when  the  complement  of 
each  sequence  is  added  a  0  (mod  4)  distribution  results. 
It  remains,  however,  that  the  complexity  of  the  generated 
sequences  are  not  determined  apriori.   The  technique  pre- 
sented is  interesting  in  its  own  right,  and  it  is  hoped 
that  the  interested  reader  may  be  able  to  apply  it  towards 
a  solution  of  the  0  (mod  4)  distribution. 

A.   RC  SEQUENCE  GENERATOR  TECHNIQUE 

Table  6.1  shows  a  listing  of  four  RC  sequences  and  their 
generators.   The  technique  used  was  to  choose  in  turn  each 
of  the  type  RC  numbers  to  have  a  functional  value  of  1. 

Note  that  the  complexity  of  S,  and  Sr  are  both  23, 
while  the  complexity  of  S,fl  and  S,2  are  both  29.   At  present 
there  is  no  way  of  knowing  which  pair  of  RC  sequences  will 
have  the  same  complexity,  except  that  it  appears  C(S.)=C(S  t-) 
iffall  four  RC  sequences  have  the  same  complexity.   A  listing 
of  the  generators  for  each  of  the  32  pairs  of  RC  sequences 
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TABLE  6.1 
RC  GENERATOR  TECHNIQUE 

C(S3)=23    C(S5)=23    C(S1Q)=29  C(S12)=29 

G(S3)       G(S5)       G(S1Q)  G(S12) 

0        111  1 

1111  1 

2        0          0          0  0 


n 


|3 

1 

0 

0 

0 

4 

1 

1 

1 

1 

IS 

o 

1 

0 

0 

6 

1 

1 

1 

1 

7 

1 

1 

1 

1 

8 

1 

1 

1 

1 

9 

1 
0 
0 

1     o 
1 

1 

0 
0 
0 
1 

1 
1 
0 
0 

1 

1 

10 

0 

11 

0 

12 

1 

13 

1 

14 

1 

1 

1 

1 

15 

1 

1 

1 

1 

indicates 

RC  number. 

S3  =  11111001010001001101110101100000 
S5  =  11111001011000001101110101000100 
S1Q=  11111001010110000011011101000100 
S12=  11111000001101110101100101000100 
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is  prepared  in  Table  6.2.   These  are  listed  in  order  by 
W(G)  and  then  by  complexity  as  far  as  the  groupings  allow. 
Keep  in  mind,  that  the  complement  of  each  RC  sequence  is  not 
included  in  the  table.   A  close  analysis  of  this  table  then 
shows  a  0  (mod  4)  grouping  for  the  number  of  RC  sequences 
of  complexity  c  and  span  n  =  5. 

B.   PROBLEMS 

For  n  =  7,  there  exist  6  type  R  numbers  (other  than  0 
§  64)  and  8  type  RC  numbers  in  the  truth  table.   This  large 
number  of  possibilities  will  create  considerable  problems 
in  extending  this  idea  to  n  =  7  and  beyond. 

Using  this  approach  a  large  scale  computer  based  analysis 
is  essential  to  any  further  investigation  of  RC  sequences 

for  n>_7 .   Considering  the  number  of  de  Bruijn  sequences  for 

57  17 

n  =  7  to  be  2   or  1.42  x  10    an  algorithm  for  generating 

17 
only  the  RC  sequences  is  crucial,  since  generating  10 

de  Bruijn  sequences  is  not  technically  feasible.   Though 

Theorem  7  establishes  the  lower  bound  of  311,040  RC  sequences 

for  n  =  7,  this  is  a  very  weak  bound  and  a  presumably  still 

conservative  estimate  would  be  on  the  order  of  10   actual 

RC  sequences  for  n  =  7. 
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TABLE  6.2 

RC  de  BRUIJN  SEQUENCE  0  (Mod  4) 
GROUPING  FOR  n=5 


G 

(S) 

W[G(S)3 

C(S) 

100  0 

1 

1 

1001 

0 

0 

1 

101 

7 

25 

0 

1 

7 

25 

1 

0" 

7 

27 

1 

ff 

7 

27 

100  0 

0 

0 

1011 

1 

IT 

0 

0 

1 

Oil 

7 

7 

27 
27 

1 

0 

7 

29 

u 

i 

7 

29 

101  0 

0 

0 

i 

0010 

0 

1 

0 

1 

1 

0 

Oil 

7 
7 
7 

31 
31 
31 

1 

ff 

7 

31 

100  IT 

1 

0 

1 

0010 

0 

1 

IT 

0 

1 

111 

7 
7 
7 

31 
31 

31 

1 

U 

7 

31 

101  0 

0 

0 

1011 

1 

0 

1 

0 

1 

Oil 

9 
9 

31 
31 

1 

ff 

9 

31 

0 

1 

9 

31 

100  0 

1 

0 

1011 

1 

0 

0 

0 

1 

111 

9 
9 

31 
31 

1 

IT 

9 

31 

0 

1 

9 

31 

111  0 

0 

0 

1111 

1 

0 

1 

0 

1 

Oil 

11 
11 

23 
23 

1 

0 

11 

29 

U" 

1 

11 

29 

101  0 

1 

T 
0 

1011 

0 

1 

0 

1 

111 

11 
11 

25 

25 

1 

(T 

11 

27 

IT 

1 

11 

27 
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APPENDIX  A 

ALGORITHM  FOR  GENERATING  LEAST  ELEMENT 
FOR  EVERY  PURE  CYCLE 


This  algorithm  is  designed  to  produce  the  least  element 
of  each  pure  cycle.   The  least  elements  are  initially  in 
decreasing  order.   A  simple  rearrangement  will  group  the 
least  elements  by  weight  in  decreasing  order  if  necessary. 
The  algorithm  to  be  presented  is  an  adaptation  of  the 
9 -operation  discussed  in  a  paper  by  Fredricksen  and  Maiorana 
[Ref.  10],  to  generate  a  lexicographic  list  of  necklaces. 

ALGORITHM:   Begin  the  0,      operation  with  the  zero  n-tuple, 
which  is  the  first  least  element  for  the  zero  pure  cycle. 

ft  Operation:   ft  (x,  •••  x  )  =  (y,  •••  v  )  =  Y 

r  v^      nl      n 

1.  Find  the  largest  subscript  j  such  that 

x.  =  0  and  x,  =  1  for  k>j  . 
J  k  J 

2.  Form  x,  x2  •••  x-i  1  where  x.  =  1 

3.  Repeat  x,  x2  •••  x-  ,1  until  n  numbers  are  produced, 
i)   If  n  =  tj  for  some  integer  t  then  let 

Y=x,  •  ••x._,lx,»««x._-,l«'*x,»»*x.,l 

ii)  If  n>tj  then  finish  with  x.,*»»x   .  , 
J  1    n"t-i 

then  let  Y=x, • • »x .  , 1* • «x, • • *x .  , lx, • • #x     . 

4.  Y  is  a  least  element  iff  n=tj 

5.  Repeat  beginning  at  step  1  until  the  final  least 
element  (1)   is  reached. 
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TABLE  A.l 

8  ALGORITHN 

[  FOR  N 

=  7 

wt 

0000000 

0 

0000001 

1 

0000010 

* 

0000011 

2 

0000100 

* 

0000101 

2 

0000110 

* 

0000111 

3 

0001000 

* 

0001001 

2 

0001010 

* 

0001011 

3 

0001100 

* 

0001101 

3 

0001110 

* 

0001111 

4 

0010010 

* 

0010011 

3 

0010100 

* 

0010101 

3 

0010110 

* 

0010111 

4 

0011001 

* 

0011010 

* 

0011011 

4 

0011100 

* 

0011101 

4 

0011110 

* 

0011111 

5 

0101010 

* 

0101011 

4 

0101101 

* 

0101110 

* 

0101111 

5 

0110110 

* 

0110111 

5 

0111011 

* 

0111101 

* 

0111110 

* 

0111111 

6 

1111111 

7 

*   Unacceptable 
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Table  A.l  illustrates  the  algorithm  for  n-7.  Arranging 
in  order  by  weight,  the  least  element  for  each  pure  cycle  is 
tabulated  in  Table  A. 2. 


TABLE  A. 2 
PURE  CYCLE  LEAST  ELEMENTS  (n=7) 


0 

0000000 

0 

1 

0000001 

1 

2 

0000011 

3 

0000101 

5 

0001001 

9 

3 

0000111 

7 

0001011 

11 

0001101 

13 

0010011 

19 

0010101 

21 

4 

0001111 

15 

0010111 

23 

0011011 

27 

0011101 

29 

0101011 

43 

5 

0011111 

31 

0101111 

47 

0110111 

55 

6 

0111111 

63 

7 

1111111 

127 
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